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Abstract: The cosmological dynamics of a brane world scenario where the bulk action is taken 
as a generic function of the Ricci scalar is considered in a framework where the use of the Z2 
symmetry and Israel junction conditions are relaxed. The corresponding cosmological solutions for 
some specific forms of f(TZ) are obtained and shown to be in the form of exponential as well as 
power law for a vacuum brane space-time. It is shown that the existence of matter dominated epoch 
for a bulk action in the form of a power law for 1Z can only be obtained in the presence of ordinary 
matter. Using phase space analysis, we show that the universe must start from an unstable matter 
dominated epoch and eventually falls into a stable accelerated expanding phase. 
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1. Introduction 

In 1999, Randall and Sundrum (RS) [Q proposed a model in which our AD universe, the brane, 
is embedded in an AdSs bulk space. All gauge fields are confined to the brane while gravity can 
propagate into the extra dimension and thus into the bulk space. As is now well-known, such a 
proposal opened a new window through which a fresh look at the universe has become a possibility. 
One outcome of the RS model was that it offered a rational explanation to the question of hierarchy; 
the enormous disparity between the fundamental forces. Since its appearance, numerous scenarios 
based on the RS setup have been proposed, each dealing with some particular aspects of the 
structure and evolution of the universe. An immediate issue after the introduction of the RS 
models was the question of what the form of the field equations on the brane is, knowing their 
form in the bulk space. An answer came along in the elegant work of Shirumizu, Maeda and Sasaki 
(SMS) |0| where they showed how to, as it were, project the field equations in the bulk onto the 
brane in a covariant manner. The subsequent avalanche of research works have been somewhat 
overwhelming. Numerous papers have appeared and tried to remedy the shortcomings of, generally 
speaking, the standard model of gravity, including the late-time acceleration of the universe, the 
galaxy rotation curves, the virial mass discrepancy, etc For a review on the brane word models 
see §. 

In the SMS method, the bulk geometry is projected onto the brane using the Gauss-Coddazi 
equations and Israel junction conditions ||. By means of the Gauss-Coddazi equations, the intrinsic 
curvature of the bulk is related to the intrinsic and extrinsic curvatures of the brane. The Israel 
junction conditions then relate the extrinsic curvature of the brane to the energy-momentum content 
of the model. In this method, the brane equations of motion reflect the global geometry of the bulk 
space-time through the electric part of the Weyl tensor. However, there are situations where the 
use of the Israel junction conditions become somewhat restrictive. In cases where more than one 
extra dimension is involved their use is not well understood. Also, there are certain matters that 
arc not compatible with these junction conditions [Bj. In the latter, one cannot relate the extrinsic 
curvature to the matter content of the brane. Such considerations motivated the idea presented in 
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and further developed in [gj where the extrinsic curvature from the Gauss-Codazzi equations is 
calculated geometrically and the assumption of Z2 symmetry is relaxed. However, the effects of the 
global structure of the bulk space is retained by the field equations. This model has been largely 
investigated in the literature ||. 

Somewhat parallel to the brane-world development, f(lZ) theories of gravity appeared on the 
scene to deal with some of the shortcomings of the cosmological implications of the standard general 
relativity. In such theories, one modifies the Einstein-Hilbert action by replacing the Ricci scalar 1Z 
with a generic function f{7Z). It is now well-known that such theories are equivalent to scalar-tensor 
theories. In fact, both the metric and Palatini formulations of f(TZ) gravity can be rewritten in 
terms of certain versions of the Brans-Dicke theory |l(| . These theories with many variations have 
been successful, for example, in describing galaxy rotation curves which is generally believed to be 
the result of the existence of dark matter. For a review on f(R) gravities see JTT[ . 

Within the context of a brane-world scenario, f(7V) theories have been studied recently in 
where the authors follow the SMS procedure to project the bulk field equations onto the brane 
hypersurface. The field equations thus obtained show that the matter content of the brane-world is 
not coupled to the bulk f(lZ) term. In this paper, we generalize the work [jl2| to a theory where the 
use of Z 2 symmetry and the Israel junction condition is relaxed. As it turns out, since the junction 
conditions are not used, the field equations on the brane contain a new tensor Q^ v which is made 
out of the extrinsic curvature and should be calculated through the Gauss-Codazzi equations. The 
effect of the non-trivial bulk action is reflected in a new tensor which is made of the bulk f(lZ) 
function. One must only project the scalar f(7V) in an obvious manner onto the brane in order to 
calculate this tensor. The resulting field equations show that the matter content of the brane is 
non-minimally coupled to the bulk f(TZ). Therefore, in contrast to the work presented in the 
non-trivial bulk action can influence the brane energy-momentum tensor directly. 

Recent observations provided by the supernovae legacy survey data Jl3| shows that, taking 
lode — const., the equation of the state parameter is found to be lode = —1-04 ± 0.06. This data 
shows that the dc-Sitter accelerating phase is a very good candidate for the late time acceleration 
of our universe. In this sense, one may expect that a good cosmological model should predict some 
matter dominated epoch followed by an accelerated expanding phase. In this paper we consider 
this issue for the brane world model with general bulk action using the phase space analysis. 
The method of phase space analysis for cosmological models is very well known and is used for 
many alternative theories |L4[ [l?]]. As it turns out, the model can explain both the matter 

dominated and accelerated expanding phases, resulting in a universe which starts from an unstable 
matter dominated phase and eventually falls into the accelerated expanding phase. 

The paper is organized as follows: In the next section we construct the model from a general 
bulk action in d dimensions, relaxing the use of Israel junction condition and the Z2 symmetry. 
The model is then restricted to a 5D bulk and cosmological solutions are obtained in section ||. 
In section ^ a dynamical system analysis is presented which shows that the model can in principle 
explain a well founded cosmological evolution. Conclusions and final remarks are drawn in the last 
section. 

2. The model 

We consider a pseudo-Riemannian manifold V m in which the background manifold V4 is isometrically 
embedded by the map y : V4 —> V m such that 

Q A By%y B v = g^ GA B y A ^ = o, q A b^K = ^ = ±1. (2.1) 

where Qab{<)hu) is the bulk (brane) metric ,{y A }({x^}) are bulk (brane) coordinates and Af^ are 
(m — 4) unit vectors, orthogonal to the brane. In order to investigate the effects of the general bulk 
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geometry, we assume a non-trivial bulk action of the form 

I = i / d m xV=Q[f(n) - 2A] + I M , (2.2) 

ZK m J 

where 

/ M = [d m x^G£M, (2.3) 



and £jvj is the matter Lagrangian and 1Z is the bulk Ricci scalar. After variation of I with respect 
to the bulk metric, Gabi we obtain 

f(n)n AB - \g AB f{n) + g AB uf{n) - v A v B f(n) + \g AB = «? m T AB . (2.4) 

where prime represents derivative with respect to the argument. Rearranging the above equation, 
we obtain the effective Einstein field equations in the bulk 

G AB = Ti AB — -jR~Q AB = S A b, (2.5) 



where 
S AB 



f'(K) 



K 2 m T AB - Ag AB - (^Kf(K) - + uf{n))g AB + V A V B f'(K) 



(2.6) 



As was mentioned in the Introduction, we use the method introduced in J?J to project the bulk 
geometry onto the brane. In this sense, perturbations of t4 in a small neighborhood of the brane 
along a generic transverse direction £ = £ a AT a (a = 1, 2, m — 4) orthogonal to the brane are 
considered as follows 

z A (x^n^y A + (^y) A , (2.7) 

where C represents the Lie derivative and £ a denotes a small parameter along M A that parameterizes 
the extra noncompact dimensions. The presence of tangent components of the vector £ along the 
submanifold V4 can cause some difficulties because it can induce some undesirable coordinate gauges, 
but, as is shown in the theory of geometric perturbations, it is possible to choose this vector to be 
orthogonal to the background . 

Let us now consider the perturbation of the embedding map along the orthogonal extra dimen- 
sion M a , giving local coordinates of the perturbed brane as 

z*(x»,e) = y%+eK»( xV )- ( 2 - 8 ) 

Using the above assumptions, the embedding equations of the perturbed geometry are given by 

q v * = QabZ* v z» v , g^ = g AB z A ^u B a , g AB M A M? = 9ah . (2.9) 

Now, use of the embedding equations and the relation 

g AB = yp B v r v + uM B 9 a \ (2.10) 

allow us to write the metric of the bulk space in the following form 



v A vb g a b ' 



- 3 - 



where 

g^u = fa - 2£ a AVa + C^g^K^aK^b, (2.12) 

is the metric of the perturbed brane and is the extrinsic curvature of the original brane defined 
as 

K^a = — Qa B %;v ■ (2-13) 

We also use the notation = ^A^d where 

A, cd = g AB M A d .^ c B , (2.14) 

represent the twisting vector fields. The presence of gauge fields A^ a tilts the embedded family of 
sub-manifolds with respect to the normal vector J\f A . According to our construction, the original 



brane is orthogonal to the normal vectors Af A . However, from equation (2.9), it can be seen that 
this is not true for the deformed geometry. Let us introduce 

K = zi - g ab M A A» b . (2.15) 

One can easily verify that for the set {X A ,Af A } we have the following projection relations 

GabX*X* = g^, GabX a JV b =0, GabM a M£ = g ab = ±1. (2.16) 

These define a new family of embedded manifolds whose members are always orthogonal to M A . 
This new embedding of the local coordinates can be suitably used for obtaining induced Einstein 
field equations on the brane. The extrinsic curvature of a perturbed brane in the coordinates 
{X A ,N A } becomes 

K^a = -QabX a N b 



— Kfiisa £ -Kfi-ya-Kvpbg^^ 

2 d£, a ' 



(2.17) 



which is the generalized York relation. In the basis {X A , Af A }, the Gauss, Codazzi and Ricci 
equations are given by 

Ra/3-fS — 2g ah K a ^ a Ks}f3 b + TZ A bcdX a X b X^X^ , (2-18) 
[jc-,8] — 2g a A[ iac Ks] a b + TZabcdX^N c X n X s , (2-19) 

and 

2^4[ / ia&;i,] = -2g m " A^ ima A v -\ nh - g ap K^ aa K^ pb - TZabcdX^X® N b N b , (2.20) 
where TZabcd and Rap^s are the Riemann tensors for the bulk and the brane respectively. After 



contracting equation (2.18) and using the relation 

G AB = X A X B g^ + Af A JV B g ab , (2.21) 

one obtains 

= (K» ac K« c - K C K^) + K AB X A X B - g ab ll ABCD N A X B Af b c X° . (2.22) 
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So, the Einstein tensor on the brane is in the following form 

Gfiv = G AB X*X B V + Q^ v + g ab K AB M*N?g^ - g ab n ABCD Af*X B Af b c ' X D V - ±Cg,y, (2.23) 



where 

O .... = -n ab ( K~l K....1. - K.. K....A A- 



Q»u = -g a \Kl a K lvh - K a K^ b ) + l(K aPa K al3a - K a K a )g^, (2.24) 



and 

C = g ab g cd n ABC D^K B K C ^F ■ (2-25) 
Using the decomposition of the Riemann tensor into the Ricci tensor, Ricci scalar and Weyl tensor 

2 2 

H-abcd = C AB cd + t 7a(Qb[d'R>c]a - Ga[dT^c]b) + 7 7w 'a\dQc]b) i (2.26) 

(to — 2) 1 ' 1 J [m— l)(m — 2) 

one obtains the generalized Einstein equations on the brane as 



G^u — G AB X^X B — £ M1/ + Q^y — -Cg^y + _ ^% g ab TZ AB J*/^J\f^ 'g^y 



where 



(to — 4) A B ( m ~ 4) 

7 t^T^abX X + - — ^Tlg^u, (2.27) 

(m — 2) ,M ' (to— 1)(to — 2) 



g^CABCDAffxZNfX?, (2.28) 



is the electric part of the Weyl tensor C AB c d ■ 

In order to obtain the field equations corresponding to action ( |2.2| ) , we follow Dvali and Shifman 
fjjll which have proposed a mechanism to localize the standard model gauge fields to the brane. 
Using this idea, we may decompose the components of the energy-momentum tensor of the matter 
as 

«4V, = -^jT/u,, T ab = 0, T M „=Q. (2.29) 
to — 2 

However, the geometric part of the generalized energy-momentum tensor includes the bulk-bulk 
and the bulk-brane components, and one must take this into account when writing the Einstein 
field equations. We thus define the components of the total energy-momentum tensor in the basis 

S^ = S ABX ^X% S IM = S A Bxfrf?> S ab = S AB N*N b B . (2.30) 
After contracting equation ( ]2.5[ ), one obtains 

1Z = 2 —S, (2.31) 

to — 2 

where S = S A A is the trace of the generalized bulk energy-momentum tensor. The 7Z AB is then 
given by 

TZ AB = t:GabS + Sab- (2.32) 

to — 2 



Substituting 1Z AB and 1Z in equation ( 2.27 ) and using equation ( 2.3C ), we find 
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Guy — Qui/ c\^9ttV ~l~ ^Sril/ I . 

2 m—2 m—2 



'""Vg^- ^l^ a (2-33) 



(m - 2)(m - 1) 

Finally, using equation ( [2.29 ), we obtain the effective Einstein equations on the brane as 



LIU Q uv ~t~ & liv 



(2.34) 



where Q M „ is defined by equation (2.24) and 



n 



w 9ltu 



2 v A v B f\n) vAvB 



^-"^ TO „2 f'(U) 



w 9fW 



2 VjNufQt) 
m-2 f'(K) ' 



(2.35) 



is a new tensor that reflects the effects of the non-trivial bulk action on the brane. In the second 
line of the above equation we have assumed X^ = 6^, so that Vp is the /i = 0,1,2,3 components 
of the m dimensional covariant derivative. We also define 



W = 



m 2 — 5m + 10 



K 



m 2 -7m + 14\ /2A-/ 



2m 2 -6m + 47 '" ' V 2(m - 2) /V /' 

m 2 -8m + 17 ^ m D/' x m-3V A V 5 f V) 



+ 



(2.36) 



m-2 y /' m-2 /' 
These tensors must be evaluated on the brane. We also note that in the limit f(7l) = 71 the tensor 



ILj„ can be converted to the trace of brane energy-momentum tensor, and the field equation (2.34) 
reduces to [Q. 

Let us focus our attention on the five dimensional bulk. In the case of a co-dimension one 
brane wc have C = due to the symmetries of the Riemann tensor. We also consider the case where 



Xj t =6$. In this case, equation (|2.21|) reduces to 



Gab = 9ab +NaNb, 



(2.37) 



and the field equations ( 2.34 ) become 



where k 2 = g/c§, and 



n MI / = Wguv - L 



(2.38) 



(2.39) 



w = T2 n + L + W 



5 Df - / - 2A 



(2.40) 



with 



_2 V M V„/' 



(2.41) 



? M „ = - X^X?) + — {K a pK a ^ - K 2 )g^, 



(2.42) 
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where K = g^ u K^ v . In the case of a constant curvature bulk, the tensor Q^ v is an independently 
conserved quantity as can be seen easily from the Codazzi equation ( 2.19 ). However, we are inter- 
ested in a general bulk geometry, so we consider Q^ v as a general tensor which reflects the bulk 
effects on the brane equations of motion. By defining the new tensor 

and using the Bianchi identity, we find that M^„ is conserved. It is therefore possible to consider 
M^ v as a new effective energy-momentum tensor made of the standard conserved matter energy- 
momentum tensor and terms reflecting the effects of the extra dimension and also the non-trivial 
bulk action. The field equation on the brane can now be written as 



(2.44) 



3. Cosmological Solutions 

To study the time evolution of the universe, we take the bulk metric as 

dr 2 



ds 2 B 



-dt 2 + a(t,wY 



1 — kr 2 



r 2 dfl 2 



dw 1 



(3.1) 



Taking M A = 8£ and using the equation ( |2.37| ) the brane metric reduces to the usual FRW form 

dr 2 



ds z 



-dt 2 + a(ty 



1 - kr 



r 2 dn 2 



With this choice, the electric part of the Weyl tensor becomes 



to ~ 2a 2 



= -3£, i= 1,2,3, 



(3.2) 



(3.3) 



—ad + a + k + ad — d 

where dot represents derivative with respect to t and a hat over an arbitrary function f(r,w) is 
defined as f(r) = §^ \ w =o- We can also calculate the extrinsic curvature and Q^ u by use of 
equations ( [2.17 ) and ( 2.42 ) respectively, with the result 



- diag(0,l,l,l), 
a; 



and 



Q fl u = ~[ -) diag(3,l,l,l). 
By assuming that the matter on the brane has the perfect fluid form 

rt = diag(-p,p,p,p) , 
the conserved tensor M^ v can be written as 

Ml = dia,g(-p A ,p A ,p A ,p A ) , 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



which, using equations ( 2.39 ), ( 2.43 ), ( |3.5| ) and (|3.6|), result in an effective energy-density and 
pressure on the brane 

2 



PA 



PA 



+ £° + W-L Q a , 
-\£%-W + L\. 



(3.8) 
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Now, with the aid of equation ( 2.41 ), L M „ becomes 



,o __2/' 



1,2,3. 



Using equation (2.44), the Friedmann and the Raychaudhuri equations become 



H" = - 



^ P + 3(-) +£0 +w + _/.__ 



a 
a 



| (p + 3p) + 2 (£° - _ jLL ( 3 /'a - 3/'a - /'a) 



(3.9) 



(3.10) 



/' vr " ^ u ' 3 a/' 
3.1 The case /(ft) = f Q e aK 

Let us consider first the case f(TZ) = /oe aK where Jo and a are some constants. Assuming = 
and A = 0, and using the trace of the field equation ( 2.44 ) 



R-K + Q 



a 



where R is the brane Ricci scalar, L^u and W take the following forms 



(3.11) 



(3.12) 



W = 



12 



K + 4(R + Q) 



With these assumptions, equation (3.13) reduces to 



W=-- 



a \a 




1 

12c? 



(3.13) 



(3.14) 



where we have assumed k = 0. The solutions of equations (3.10) are of an inflationary form, given 
by 



a{t) = c 2 e 
a(t) = /3a(t) 



Ht) = -(p 2 + — )a(t) ) 



(3.15) 



where c\ and c 2 are constants which we choose to be positive and, (3 = /3(a) is an arbitrary constant 
depending only on a. As one can see from the above solutions, derivatives of the 51? scale factor 
with respect to the extra dimension w is proportional to the 4D scale factor. 

We should note that in our derivation of the field equations for the brane, we did not use the 
Codazzi equation. However, the Codazzi equation must be satisfied for all smooth manifolds. One 
can easily check that the above cosmological solutions satisfie the Codazzi equation. 

3.2 The case f(K) = f (K - K ) a 

Now, we consider a power law form f(lZ) = fo(lZ — ~R-o) a , where /o, 7?-o and a are some constants. 
Using the trace of the field equations ( 2.44 ), one obtains 



Q-(l--)K 



no 

a 



(3.16) 
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Next, using ( |l40| ) and (HH), one finds 



W 



1 

12 



(3.17) 



(a - 1) (a - 2)VJlV v 1l + (K- K Q )V^„K 



(3.18) 



For this specific form of f{TV), we again consider the solutions of the equations ( |3.10| ) in the absence 
of ordinary matter and cosmological constant. For general a one has a solution 



a(t) = at 1 / 2 , a(t) = c 2 a(t), d(t) = - ( cf. + -K ) a(t) 



6 



(3.19) 



where c\ and C2 are the constants of integration. This form of the scale factor is similar to the 
radiation dominated epoch. 

We have also a power-law solution 



a(t) = Cif 



6(t) 



-a(t), a(t) 



1 



U a(t), 



(3.20) 



where 



2a 3 - 9a 2 + 4a + 6 
a 2 + 7a - 12 ' 



(3.21) 



and ci and C2 are the constants of integration. This solution shows that the above form for the bulk 
action cannot admit a matter dominated behavior. However, this can be traced to our assumption 
of the space-time being empty. In the next section we will see that adding ordinary matter to the 
system will produce matter dominated solutions. We have also an inflationary solution as follows 



<i(l) = r l <" 1 . ,-,(/ ) = <- 2 ,i[l i. »/(/) = ( 2n 2 -4~^o j «!/!• 



(3.22) 



where n, c\ and c-i are some constants. The solutions above also satisfy the Codazzi equation and 
are therefore cosmological solutions of the model. 

All the solutions above have the property that derivatives of the scale factor with respect to 
the extra dimension are proportional to the scale factor itself. In the next section we will use this 
property to build an autonomous system of equation for the theory. 



4. Cosmological Dynamics 

In this section, we investigate the cosmological dynamics of the model introduced in the previous 
section. Substitution of expressions for the electric part of the Weyl tensor and quantity W in 



equations (3.10) results in the Friedmann and Raychaudhuri equations for a general f(lZ) 



-H (J +- + l! + 2— -2---( Prad + Pm )=0, (4.1) 



and 



where we have defined F(7Z) = f'(TZ). With definition of the bulk Ricci scalar 



K = 6 



2H 2 + H 



one can write the Friedmann equation as 



1 = 



n 



(A_\_l_f i_ 



Fa 



12H 2 \HaJ 8 FH 2 HF H 2 Fa 2FH 2 



(Prad + Pm) ■ 



(4.3) 



(4.4) 



In order to study the dynamical evolution of the system, we define the following dimensionless 
quantities fl4| 



and 



1 



1 

12' 



■r-2 



1 / 



H 2 F' 



x 3 



F Fa 

£4 



HF 



H 2 Fa' 



K Prad , . r\ 

X5= 2IPF> (4 ' 5) 



a2 , i 



f^Pm 

2H 2 F 



^rad — 



Qdark — ? i 
i=l 



(4.6) 



with physical constraints represented by f2 m ,Sl rQ( 2 > 0. As was mentioned before, the Codazzi 
equation must be satisfied for any solution to the field equations. In our case, it is reduced to 



3° + J=0. 
a b 



(4.7) 



We now assume that the derivatives of the scale factor with respect to the extra dimension are 
proportional to the scale factor itself, that is 



a - = b, 



a 

- = c, 
a 



(4.8) 



where b and c are constants which may depend on the detailed functionality of f(JZ). This, for 
example, occurs in the cosmological solution obtained in the previous section. We also define the 
dimensionless quantity 



x 6 = 

The resulting dynamical equations become 

dxi 



K + 6b (b + c) 
K + 12b 2 ' 



A X l X 3 , 2 

dA 3m 



(4.9) 



(4.10) 



dx 2 x x x 3 
_=4.* 2 + W6 + — 



+ x 2 x 3 , 



(4.11) 



-j^ = -1 + (2/3 +l)x 1 + ^ftxxxa - 3x 2 - 2x 4 + x 5 + x\ + ^ Xl x 3 x 6 , 



(4.12) 



= - (2/3 + 6) xi - i (3 + /3) xix 6 + 4.T4 + xiX^Xq + £3X4, 



(4.13) 
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dx 5 
dA 



XiX 5 X 6 + X 3 X 5 , 



and 



where A — In a and we have defined 



We also define the quantity 



dx 6 1 . „. 
dZ = 3^ X3(:C6 + 4) ' 



P = 



3c 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



which depends on the bulk Ricci scalar. From equation (4.9) we see that the Ricci scalar can be 
represented as a function of the quantity xq , so one may consider the quantity m as a function of 
x 6 . 

Let us define the effective equation of state ptot = ^effPtot, where 



2_H_ 

3H 2 



xix e ) . 



(4.18) 



With these definitions, w e ff = gives the behavior a oc frs which is the standard behavior of the 
matter dominated epoch. Also ^ e ff = — 1 corresponds to the de Sitter epoch. 

In order to find the dynamical behavior of the above cosmological model, one must obtain the 
critical points of the dynamical system ( 4. ld| )- ( 4.15 ) followed by expanding the system near the 
critical point. The resulting equations can be written as 



X' = EX. 



(4.19) 



where X is a column vector made out of Xj's, and the prime represents derivative with respect to 
parameter A = In a while 5 is a matrix obtained from linearizing the system ( 4. 10| )- (4.15). The 
eigenvalues of this matrix at each critical point determine the stability of that point. If all the real 
parts of the eigenvalues are negative, then the point is stable. The appearance of positive eigenvalues 
makes the point a saddle point. If, however, all real parts of the eigenvalues are positive, we have 
an unstable point. In the evolution of the universe, one needs a saddle or an unstable point which 
would correspond to the matter epoch. Because of the instability of the this phase, it can fall into 
a stable point which would then correspond to an accelerating phase. 



The critical points of the dynamical system (4. 10)- (4. 15) are discussed in what follows: 



4.1 Radiation epoch 



Pi : (xi,x 2 ,x 3 ,x i ,x 5 ,x e ) = (0,0,0,0, l,x) 

"eff = 5, 



a m = o, n rad = i, n dark = o, (4.20) 

where x is arbitrary. This is the standard radiation dominated epoch, with a oc t^, using equation 



(4.18). The eigenvalues of this point are 

0, 1, -1, 4, 4, 4, 
making this point is a saddle point as discussed above. 



(4.21) 
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4.2 Kinetic epoch 

We have three critical points which correspond to the kinetic epoch of the universe as follows 

P 2 : (xi, x 2 , x 3 , x 4 , x 5 , x 6 ) = (0, 5, -4, 0, 0, -4) , 

W e// = = 0, VLrad = 0, ^dark = 1- (4.22) 

P3 : (xi,X2,x 3 ,X4,X5,x 6 ) = (0,0,1,0,0,-4), 

W e// = = 0, VL r ad = 0, ^darfe = 1, (4.23) 



and 



Pi : (xi,X2,X3,X4,x 5 ,x 6 ) = (0,0,-1,0,0,-4), 

W e// = £l m = 2, O ra d = 0, fldark = — 1. (4.24) 

The eigenvalues to these points are 

P 2: „,_ 5 ,_ 4 ,_ 3 ™,_J_, ,, 25) 

P a :l,2,l?^-J,5,5,-L, (426) 
3m 3m 

P 4 :-2,-l,l^±l,3,3,--L. (4.27) 
3m 3m 

As can be seen, point P4 is a saddle point. The point P 2 is either stable or saddle, and the point 
P3 is either saddle or unstable depending on the sign of in. 

4.3 de Sitter epoch 

This epoch includes four critical points corresponding to the accelerating phase of the universe. 
The one which can be used to describe the late time acceleration of the universe is the stable point. 
However, unstable accelerated fixed points can be used to describe the inflation epoch to ensure the 
end of the inflationary accelerating phase. The critical points corresponding to this epoch are 

/ 3m -9m 12m 144m 2 + 21m -2 \ 

P 5 : Xi,X 2 ,X3,X4,X 5 ,X 6 = — — — , — —2 , — - — ,0, — —2 ,-4 

^1 + 3™ 2(1 + 3m) l + 3m 2 (1 + 3m) / 

l-9m _6m(l + 6m) 144m 2 + 21m -2 _ 3m (30 m + 7) 

3 + 9m (l + 3m) 2 2 (1 + 3m) 2 (1 + 3m) 

(4.28) 

For m < — -i this point corresponds to the phantom epoch, and it is also compatible with the 
condition Q m > 0, but it gives fi ra d < 0. The corresponding eigenvalues are then given by 

4 ax + 4m7 a x - 8mj ± i\ f Za 2 (4 29) 

' l + 3m' 27(1 + 3m)' 4(1 + 3m) 7 ' [ ' ' 

where we have defined 

2 

ai, 2 = (-2240m 3 - 1572m 2 - 156m + 16 + 2s/m} 3 ± (208m 2 + 14m - 6) , (4.30) 

and 

7 = (-2240m 3 - 1572m 2 - 156m + 16 + 2\/m) 3 , (4.31) 
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with 



M = -995328m 6 + 1306368m 5 + 956628m 4 + 130218m 3 - 11226m 2 - 1626m + 118. (4.32) 

The critical point P5 is thus a saddle point and together with the result fl ra d < cannot serve as 
a final accelerating phase for the universe. The next point is 

P 6 : (xi, x 2 ,x 3 , Xi,x 5 ,x 6 ) = (1,0,0,0,0,-4), 

LO e ff — —1, Q m = 0, Qrad = 0, Qdark = 1, (4.33) 

which corresponds to the standard de Sitter epoch. The eigenvalues corresponding to this point are 

1 /3 

0, 0, -4, Si -2, --6i-2±i^-S 2 , (4.34) 

where 

ft ,,_JL ± 2i!»^i, 7 _ m l UJE-J)', (4.35) 
om 37 \ V m I 



and 



M = 2662 - 1548m + 38016m 2 - 27648m 3 . (4.36) 



So for —4 < Si < 2 we have a stable de Sitter epoch which, when written in terms of m, leads to 
m < —A, 0<m<iiorm> 1.38. The point P$ corresponds to an accelerating epoch much in the 
same way as the cosmological constant. As was mentioned in the Introduction, the observational 
data predict that the accelerating behavior of the late time universe is very close to that of the 
cosmological constant. Our model is therefore compatible with such predictions. Finally we have 
two other critical points which can be written as 

Pr.s ■ (xi,x 2 ,x 3 ,X4,x 5 ,x 6 ) = (-^-x(l + 3m),x,4m[2x(l + 3m) + 3] ,0, 0, —4^ , (4.37) 



J eff 



-x(l + 3m) + 1 



— 1 ^darki ^rad — 0; ^dark — ^ ^ 



x(l + 4m) + 2 



where 



1 11 - 42 m - 720 m 2 - 1728 m 3 ± V20736 i 
64 



m (1 + 12 m + 45 m 2 + 54 m 3 ) 



(4.38) 



15552 m 3 - 9468 m 2 - 540 m + 121 



(4.39) 



i < m < 0. In 

6 



The point P-j is interesting since it represents a phantom epoch for m < — -| c 
figure |l| we have shown the behavior of (jJ e ff as a function of m. However, we must impose the 
physical condition Q m > to this point which leads to —0.21 < m < 0. Figure |^ shows four 
different eigenvalues of the critical point P7 as a function of m. As we can see from the figure, 
only in the case — jjr < m < all eigenvalues are negative. So the critical point P^ is stable if 
— 4 < m < 0. This is interesting since such a range is also physical and can be used for the stable 



accelerated expansion phase of the universe, 
phase for |m| ^> 1 



The critical point P% has an accelerating de Sitter 
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Figure 3: Plot of w e // as a function of m for the critical point Ps- The asymptotic value for w e // in the 
limit m — > oo represents the de Sitter phase with w e // = — X. The matter dominated phase is represented 
by m — —0.36 as shown in the figure. 



Uaf/ = -1, il m = (4.40) 



The eigenvalues for the corresponding point in this limit are 



4m ' 3m ' 6m ' 4m ' 4m ' ^ 
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Figure 4: Stream plot of the point Ps and Pe which correspond to the unstable matter dominated phase 
and stable accelerating phase respectively. 



The point Pg, is thus stable provided m < 0. In figure g we have plotted w e // as a function of m for 
the point Pg. As can be seen, u e ff reaches its minimum value at minus infinity which corresponds 
to the de Sitter accelerating phase. Figure || shows another interesting feature of the critical point 
P%. The value m = —0.36 corresponds to w e // = which has the behavior of a matter dominated 
epoch, a(t) oc ti , with 



O m = 0.14, 



^rad — 0, 



n 



dark 



= 0.86. 



At the fix point, we have from (4.17) 

f(K) =/o(ft + 26 2 ) 1 . 
The eigenvalues corresponding to this point are 

-0.3, 3, -4.6, -3.5, -0.3, 



-0.3. 



(4.42) 



(4.43) 



(4.44) 



This point is therefore a saddle point and represents the matter dominated era of the universe. 

The discussion in this section suggests that the generalized model proposed here has a dominant 
accelerated expansion behavior. This is desirable since the brane-world models are well-known for 
the production of late time accelerated behavior of the universe. However as we saw in this section, 
the effect of the non-trivial bulk action results in changes in the behavior of the scale factor and 
the dominant contribution is that of dark energy. 



5. Conclusions and final remarks 

In this work we have considered a brane-world scenario in the context of f(TZ) gravity. We showed 
that the change in the dynamics of the bulk space-time leads to changes in the dynamics of the 
brane space-time. The change to the brane equation of motion is however important, for we have the 
non-minimal coupling of the brane energy-momentum tensor to the bulk f(TV), which is similar to 
4-dimensional f(R) gravity models. This would not have been the case had we used Israel junction 
conditions to project the bulk geometry to the brane fl2fl . 
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The cosmological solution found with the exponential form of f(TZ) exhibits an inflationary 
behavior. The power law assumption for the form of f(lZ) was also considered. The main point 
of this ansatz is that if we consider the vacuum space-time one cannot obtain matter dominated 
solutions for the theory. However, this difficulty can be resolved if one adds an ordinary matter to 
the theory. On the other hand, the radiation dominated solution always exists. It is worth noting 
that all the solutions obtained in this paper have the property that the derivative of the scale factor 
with respect to the extra dimension are always proportional to the scale factor itself. Of course 
solutions with general extra dimensional derivatives may also exist. 

In order to describe the cosmological evolution of the universe, one must consider the general 
form of f{TV). This was done by studying the phase space of the theory. The interesting result was 
that there is an unstable critical point Js, which can be considered as a matter dominated phase 
of the universe. This point corresponds to a power law form for f(lZ) and can fall into a stable 
accelerating critical point Pq which leads to the value m = —0.36. Figure (Q) shows a stream plot 



of these two points. We also note that with the assumption (4.8), we obtained a new dynamical 



variable which contains derivative of f(lZ) with respect to the extra dimension. This dynamical 
variable is a new feature of the model and is responsible for the change of the dynamics of the 
system relative to what one gets in the corresponding AD scenario jl4|. The complete analysis of 
the system would be more complicated if we relaxed the assumption (L8) which reflects the full 
behavior of the non-trivial dynamics of the bulk space-time. 
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